Probing Gluonic Spin-Orbit Correlations in Photon Pair Production 
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We consider photon pair production in hadronic collisions at large mass and small transverse 
momentum of the pair, assuming that factorization in terms of transverse-momentum dependent 
parton distributions applies. The unpolarized cross section is found to have azimuthal angular 
dependencies that are generated by a gluonic version of the Boer-Mulders function. In addition, the 
single-transversely polarized cross section is sensitive to the gluon Sivers function. We present simple 
numerical estimates for the Boer-Mulders and Sivers effects in diphoton production at RHIC and 
find that the process would offer unique opportunities for exploring transverse-momentum dependent 
gluon distributions. 
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Introduction. — Hard hadronic processes with small 
transverse momentum qt of an observed final-state sys- 
tem have attracted a lot of interest because of their sen- 
sitivity to intrinsic parton transverse momenta. Such 
processes may hence offer detailed insights into the par- 
tonic substructure of hadrons, in terms of transverse- 
momentum dependent parton distributions (TMDs). Of 
particular interest are correlations of the parton trans- 
verse momentum with the nucleon or quark spin, which 
are expressed by the Sivers [1] and Boer-Mulders (BM) [2] 
functions. From these, one ultimately hopes to learn 
about spin-orbit correlations and orbital angular mo- 
menta of partons confined in a nucleon. So far, the main 
focus of the field has been on quark TMDs. This is due to 
the fact that quark TMDs are primarily probed in semi- 
inclusive deep- inelastic scattering (SIDIS) and the Drell- 
Yan (DY) dilepton production process, which have been 
accessible experimentally [3]. On the theoretical side, the 
relative simplicity of these two reactions has allowed to 
derive factorization theorems involving TMDs [4. 5]. 

Gluon TMDs [6] and processes sensitive to them have 
received closer attention only quite recently, at least for 
cases where nucleon or gluon polarization matter. Several 
processes for accessing the Sivers [7, 8] gluon distribution, 
or the gluonic version of the BM function (more appro- 
priately described as the TMD distribution of linearly 
polarized gluons in an unpolarized nucleon) [9, 10], have 
been proposed for high-energy hadronic collisions, in par- 
ticular, at the Relativistic Heavy Ion Collider (RHIC), or 
for ep scattering at a future Electron Ion Collider (EIC). 
There is a generic dilemma concerning the processes con- 
sidered so far: in cases where experiments can be carried 
out at today's hadron colliders, factorization is known to 
be broken for TMDs [11], or to hold at best for weighted 
asymmetries that only give information on integrated 
TMDs with certain transverse momentum weights. On 
the other hand, while transverse-momentum dependent 
factorization is expected to hold for reactions such as 



ep — > ccX or ep — jet jetX [9], realization of an EIC is 
still a decade or so away. 

In this Letter we argue that the process pp — )■ 77X can 
be used to study spin-dependent gluonic TMDs in a the- 
oretically clean process already at RHIC. In proposing 
this process, we are motivated by the following observa- 
tions. First of all, since the final state is a color singlet, 
the diphoton process is expected to share many features 
with the DY process, as far as factorization is concerned. 
Indeed, like the DY process, its lowest-order contribu- 
tion comes from qq annihilation, qq — )■ 77, which can be 
shown to give rise to the same Wilson lines as the DY sub- 
process gg — > 7* , and hence involves the same quark and 
antiquark TMDs. Second, it has been known for a long 
time that in the spin- averaged case [12] at colliders pho- 
ton pair production is in fact dominated by the process 
99 ~^ 77i that is, gluon-gluon fusion to a photon pair via 
a quark box. Even though this process is formally down 
by two powers of the strong coupling constant as with re- 
spect to qq — >■ 77, the suppression is compensated by the 
structure of the associated hard-scattering function, and 
by the size of the gluon distribution function. Hence, an 
experimental study of gluon TMDs should in principle be 
possible in this process. Finally, in order to study TMDs, 
precise measurement of the (small) transverse momen- 
tum of a final state is crucial. It seems to us that this 
should be easier to achieve for a photon pair than, for 
example, for the jet pair in the reaction ep — > jetjetX. 

Being a background to a possible Higgs boson decay 
into two photons, QCD diphoton production has received 
a lot of attention in theoretical studies, in particular, 
for the diphoton pair transverse-momentum distributions 
based on perturbative all-order resummation of Sudakov 
logarithms [13, 14[. In fact, these studies pointed out 
that the resummation formalism naturally suggests the 
presence of gluon TMDs, among them a perturbative 
spin-flip distribution akin to the gluonic BM function. 
In our present Letter, we examine the diphoton process 



entirely from the point of view of TMD factorization. 
Focusing on the gluonic Sivers and BM functions, we re- 
strict ourselves to the application of an effective tree-level 
TMD formalism in the spirit of Refs. [15, 16]. At present, 
we are not able to present a proof that TMD factoriza- 
tion indeed holds for this process. Given the color-singlet 
nature of the final state and its similarity to DY kinemat- 
ics, it appears plausible that such a factorization could be 
established if Q ^ pj- ^ g-p. where Q (qt) is the photon 
pair mass (transverse momentum) and px the transverse 
momentum of one photon. We hope that our study will 
motivate work in this direction. 

Measurements of diphoton production have been car- 
ried out at the Tevatron [17]. Detection of diphoton sig- 
nals should be well feasible in polarized pp collisions at 
RHIC [18] - statistics for the reaction will depend of 
course on the collected luminosity. Concerning the ex- 
traction of TMDs from pp — > 77X, a potential compli- 
cation arises due to the fact that photons can also be 
produced in jet fragmentation, which would very likely 
spoil TMD factorization. Such fragmentation contribu- 
tions may be strongly suppressed or even eliminated by 
using isolation cuts on the photons. We leave a more 
detailed discussion of this issue to a future publication. 

Kinematics. — We analyze the diphoton process 
h{Pa) + h{Pii) — > 7(qa) + j{qb) + X in the center-of-mass 
(cm.) frame of the incoming hadrons with momenta 
= VW2 [0,1, Ot] and Pj^ = ^3/2 [1,0, Ot], where 
we used the light-cone notation = [0 



with 
Pb)^. 

The expressions for the photon momenta are simplified 
for < Q [16] 



= (a° ± a^)/V2, ar - ia\a^), and S = (Pa 



%{b) 



^^[xfci^,a;,i^,±Vi^sin0e^] , (1) 



where the upper (lower) sign in the above expression 
refers to photon a (b), Xa/b = Q^/i'^Pa/b ■ q) with the 
photon pair momentum q = Qa + Qb, and the spatial ori- 
entations of the photons are fixed by their angles 9, (j) 
in the Collins-Soper (CS) frame [4, 16]. The Lorentz 
transformation between the cm. frame and CS-frame 
has been worked out in Ref. [16] for the (kinemati- 
cally identical) DY process hh — > i'^i^X. In Eq. (1), 
60 = (cos(/), sine/)). Additional azimuthal dependence 
may be introduced by transverse spin vectors of the 
hadrons, Sa/bT = {cosipa/b, sm(f>a/b)- The partonic 
Mandelstam variables expressed in the cm. -frame read 



qa 



^ sin^ I and 



u = — 2fcf, • Qa ~ —Q^ cos^ |, where fc^^t, are the incoming 
parton momenta. 

Photon pair production in qq annihilation. — At low- 
est order, photon pairs are produced through quark- 
antiquark annihilation, qq ^ 77. By following the same 





Figure 1: Photon pair production by gluon-gluon fusion, 
steps of the DY calculation in [16], we find for qr ^ Q 



dcr««' 



d^qdn 



dcr««- 



qT<s:Q 



d^qdn 



qT<.Q 



(e^^e^), 



(2) 

where the expression for dcr'^'^^' ' can be found in 
Ref. [16]. In Eq. (2), the overall factor 2/sin^ 9 is caused 
by the fact that the process qq — > 77 proceeds via t and 
u channels while the DY process is via the s-channel. 
As pointed out above, the diphoton production and DY 
processes share the same quark and antiquark TMDs be- 
cause both have only initial state interactions. 

Gluon TMDs and photon pair production. — Gluon 
TMDs are defined through the correlator [6] 



1 



dz d zt f^j-k-z 
(2^)" 



(3) 



x(P,5|f;,(0)W'^''[0;z]F,^,(z)|P,5) 



z+=0 



which is a diagonal hadronic matrix element of two field 
strength tensors F^^'^ between nucleon states with large 
momentum component P^ and spin vector S. The Wil- 
son line W"'' with color indices a, h in the adjoint rep- 
resentation makes the correlator gauge invariant [5, 11]. 
However, the explicit form of the Wilson line depends on 
the color structure of partonic scattering that the gluon 
TMDs are convoluted with. For the gg — > 77 subprocess 
of photon pair production, we have all gluon TMDs with 
past-pointing Wilson lines. 

The leading terms in a 1/P+ expansion of T f^u-x^i are 
given by r+*'+^ with transverse indices i,i £ {1,2}. 
For an unpolarized ([/) or a transversely polarized (T) 
hadron of mass M one has the following decompositions 
of the correlator r+*'+J into gluon TMDs [6, 19]: 



^ U 

i J, 



(x, kr) 



2 M 
8M IT 



J : 



IT ^ 



hi', 



M 



9lT 



M "IT 



his, (4) 



where = e and 



Each 



of the TMDs in (4) is a function of x and 



The un- 
polarized correlator Tjj contains two gluon TMDs, the 
unpolarized gluon distribution ff and the gluonic BM 
function h^^ . Tt is parameterized by the gluon Sivers 
function fifl, and several other gluon TMDs [6]. Among 
these, his ^^'^ gluon "pretzelosity" TMD. 
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The leading order diagrams for gg — >■ 77 are shown 
in Fig. 1. Helicity amplitudes for them have been pre- 
sented in Refs. [20, 21]. While these could be combined 
in a way suitable for projecting onto transverse exter- 
nal gluon indices, we choose to compute the diagrams 
directly, defining "semicontracted" amplitudes M±± = 
M"^^ (£j.((7a))* {£±iqb))* , with transverse gluon indices 
and contracted with photon polarization vectors 
s±{qa/b)- To perform the loop integrals, we use stan- 
dard Mellin-Barnes space methods. The semicontracted 
amplitudes A1±± are finite. We obtain 



M% = /C(/t44 + /„44+£le*^ + e;£l) ,(5) 



with s± = (l,±ii)e=F»<^, IC = 25"''a,aomE9e^, where e, 
denotes the fraction of the charge of the quark in the 
box in units of the elementary charge e and a, b are the 
gluon adjoint color indices. Furthermore, fs = L(t,u)^ 
ft = L{s,u), fu = L{s,t), with L defined as 



L{x,y) = l-|^(ln|f|-i.0(-f) 



2 (x+yP 



(6) 



Note that we have fs = -A/++ , ft = -A/+2+_, /„ = 



-M 



(1) 

H h 



in terms of the helicity amplitudes of [21]. 



Our semicontracted amplitudes and the correlator 
p+';+j of Eq. (3) can now be used to compute the cross 
section for gg 77 in the TMD formalism: 

-jj-^ ^ nj SkaTd^hT S'-^HkaT + hx " Qt) X (7) 

Al ,A2 

where n = {128{2TT)'^XaXi,S'^)-'^ , and where we sum over 
the photon helicities. 

Using the decomposition in Eq. (4) we derive from 
Eq. (7) the following result for the unpolarized and sin- 
gle transverse spin polarized cross sections in terms of 
the CS-frame angles: 



= a^,^[:F,{e)c[ff ff] + :F2{e)c[w,hi^ hi^ 

+ cos(20){ J-3(0) (C [lui h^^ /f ] + C [w2 f! hi^^ ) } 



TU _ og 



(8) 



(9) 



where ctq^ — 2IC^'H, and where the ellipses denote ad- 
ditional terms that vanish upon 0-integration. We have 



defined ^,(9) = + [/^p + |/„|2 + 5, ^^^e) = 2(/,, - 1), 
= + Wu + ft] - 1, MO) - fuf; + ftf* + 2, and 

C[wf,f2] = J cPkaTd^ktrS^^HkaT + hT-qr) 

X w{kaT,hT)fl{Xa,PaT)h{xb,klj,). (10) 

Defining {ah)± = {aibi±a2b2)/2AP and [ab]± = (ai62± 
02^1 )/2M^, the following weights appear in (8) and (9): 

Wl = -2{kaTkaT)- , W2 = -2{kbTkbT) ^ , 
W-i = JjkaT,! , W4 = (kaTkbr)- " [kaTkbT]^ , 

W5 = [karkbr]- - [kaTkbT)\ i 

= 2l7 {{kbTkbT) + kaT,2 " '2{kaTkbT) + kbT,2) , 

W7 = -Jj{kaTkbT)+[kaTkbT]-kaTS ■ (H) 

We stress that the angular structure of the unpolarized 
cross section shown in (8) is identical to that found in the 
context of collinear factorization for perturbative soft- 
gluon radiation from the LO process gg —J- 77 [13, 14]. 
This may hint at a possible matching of the TMD 
and collinear formalisms in the intermediate qt region 
Aqcd ^ Qt ^ Q- We also note that weighted cross sec- 
tions of the form (F) = / (P qr dcj) F {qr , (j)) {da/d^qdVi) 
may help in disentangling the various terms in (8) and 
(9). For instance, (5;^ cos(4(/))) oc /if ^^^^(xq) /^^^^^^(xh), 
and (g^ cos(2(/))) cx h\^'''^^^ [xa] ff'^^ {xb), with kr mo- 
ments of ff and h^^ . 

IV. Numerical estimates. — In order to estimate the 
size of the various contributions to (8), (9), we use a Gaus- 
sian model for the TMDs as frequently chosen for the 
analysis of SIDIS or DY data ]3, 22]. For the unpolar- 
ized quark and gluon TMDs fl'^ we make the ansatz 



/i (2;, ^t) 



^c" 7 , (12) 



with widths (3 and 7 for which we assume /? = 7 = 
0.5 GeV^ at RHIC, and with the /cT-integrated parton 
distributions of ]23]. Very little is known about the 
other quark and gluon TMDs at RHIC energies. Model- 
independent positivity bounds for them were derived in 
Refs. ]6, 24]. To estimate the maximally possible effects 
in the diphoton process we assume saturation of these 
positivity bounds for both quarks and gluons. For the 
gluon Sivers function this gives approximately 



M rg 



(13) 



Similarly the positivity bounds lead to the following ap- 
proximations for the other TMDs: ''I ~ {2M'^)/k'^ /f , 
\hi'i\^ \ft^\ M/krfl M/krf? (with hf = 

h{j, + /fc2/(2M2)/if/), and \h^^\ ~ {2M^)/k^ff. 

Ill Fig. 2 we present our numerical estimates from our 
Gaussian ansatz. In generating those curves, we required 
each photon to have a transverse momentum of at least 
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Figure 2: Pair rapidity (y) dependence of the various terms 
in the cross sections in Eqs. (8) and (9), for the unpolarized 
[top] and single transversely polarized [bottom] cases, in pp 
collisions at ^ = 500 GeV. For the spin-dependent cross 
section we show the absolute value since the sign of the TMDs 
is not fixed by the positivity bounds. For comparison we also 
show predictions for the unpolarized Drell-Yan process in the 
upper panel, without any cuts on lepton transverse momenta. 



1 GeV, and we integrated over 4 < < 30 GeV^, 
< QT < 1 GeV, and the CS angles with appropri- 
ate azimuthal weightings. For the unpolarized cross sec- 
tion (upper panel), the gg — >■ 77 channel dominates at 
midrapidity while the qq — >■ 77 channel is more impor- 
tant at forward/backward rapidity {\y\ > 2) of the pho- 
ton pair. The contribution by the gluon BM effect to 
the 0- iridependent cross section turns out to be rather 
small. On the other hand, the cos2(/) and cos 40 contri- 
butions induced by gluons could be at the percent level 
for TMDs saturated by the positivity bounds. Realisti- 
cally, however, one may expect smaller effects depending 
on the actual size of the TMDs. In order to estimate 
the maximum size of quark Sivers contribution to the 
spin-dependent cross section, we kept only the positiv- 
ity bound saturated up-quark Sivers function since the 
up- and down-quark Sivers functions have an opposite 
sign. From the lower panel in Fig. 2, it is important to 
note that the gluon Sivers effect exceeds the quark Sivers 
effect by a factor seven or so at midrapidity, and domi- 
nates the contribution for a wide range of rapidity. That 
is, the single transverse spin asymmetry of the diphoton 
production at RHIC could offer excellent opportunities 
for exploring the gluon Sivers function. Other effects 
caused by the gluon TMDs h^j, and hfj? are negligible. 

V. Conclusion. — We have investigated photon pair 
production in hadronic collisions in the framework of 
TMD factorization. We have shown that this process 
may be suited for studying gluon TMDs at RHIC. The 
cos(40) modulation can be used to extract the gluon 
Boer-Mulders function. Even a small effect can be sig- 
nificant since this modulation is absent in the qq chan- 
nel. The cos(20) modulation ultimately gives informa- 



tion on the sign of h^^ . Such measurements may also be 
performed at the LHC where the production rate from 
gluon fusion is much larger. Another unique feature of 
the diphoton process is its sensitivity to the gluon Sivers 
function in polarized proton collisions. Measurements 
at RHIC could hence provide important clues about the 
correlation between gluon motion and hadron spin. 
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